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of elasticity using an observed solution to the forward problem. Under 
appropriate conditions this problem can be treated as a first order hyperbolic 
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the convergence of our numerical scheme and obtain error estimates. 
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1. INTRODUCTION 


A class of control and Identification problems for which the models are 
based on the equations for elastic structures are those dealing with large 
space antennas. Mathematical models of these problems are based on the 
partial differential equation 

< 1 -‘> h (e &> + h (e = f 

where u(x,y) is the vertical displacement of the antenna surface, f(x,y) is 
the distributed loading force per unit area, and e(x,y) is the distributed 
coefficient of elasticiy of the antenna surface [1], 

The identification of e using measured u and f values for the 
antenna surface is an important inverse problem. A common identification 
strategy is the "indirect" one in which one minimizes via an iterative process 
the deviation between a computed forward solution u g and the observations 
(see, for example [1]). Alternatively, e can sometimes be identified by a 
direct approach involving approximate solution of the hyperbolic equation 

(1.2) Ve*7u+ eAu = f, 

for example, by seeking the finite dimensional representation for e which 
minimizes the residuals of a difference approximation for equation (1.2). 
This is referred to as the "equation-error" method. 

A practical limitation to the direct approach for identifying e is that 
the coefficients of the hyperbolic problem (1.2) involve derivatives of the 
measured quantity u. However, when it is feasible it is simpler and cheaper 
than the indirect approach. 
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In [2] Richter presented a systematic analysis of the inverse problem 

(1.3) Ve • Vu + eAu = f, x ft IP, 

in which the coefficient e is to be determined on the basis of an observed 
(f,u) pair. He showed that the hyperbolic problem (1.3) has a unique 
solution assuming prescribed values along the inflow portion of 
3ft for f£L°°(ft), provided 

(1.4) inf [max{ |Vu(P) | ^ , Au(P)}] > 0. 

P£ft 

He also proved that if condition (1.4) holds, then e depends continuously 
on f in L°°(ft). 

In this paper, we show that if condition (1.4) holds then e depends 

continuously on f in L^(ft) for all p [ 1 ,°° ) . We then use this 

2 

continuous dependence result in L (ft) to propose a spline-based technique 
for approximating the unknown coefficient e in equation (1.3). We prove 
that our scheme converges to the actual solution e of (1.3) and obtain error 
estimates. 

In [2] Richter proposed especially favorable "test conditions" for 
observing a forward solution u to the elliptic problem for (1.3): 

(1.5) inf f > 0, u = 0 on 3ft. 

ft 

Under these conditions (1.4) will be satisfied and the hyperbolic problem 
(1.3) with the resulting (f,u) pair will require no Cauchy data for e 
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because the characteristics of e will originate at points of degeneracy 
within £2, rather than on 3£2. 

Typically, condition (1.5) holds for antenna problems. Thus, the 

numerical algorithm we propose in this paper, and which is based on a 

2 

continuous dependence result in L (£1), would be particularly suitable for 
antenna problems. In a companion paper [4], we present a multigrid algorithm 
for approximating e numerically. 


2. CONTINUOUS DEPENDENCE OF THE HYPERBOLIC PROBLEM 

Let B(u)e be the operator 

(2.1) B(u)e »V • (eVu) - f, x £2 C IP, 

where e and u are defined in a connected, bounded domain £2 C if 1 , 

2 — 1 — 

Throughout this paper we shall assume that u€C (£2) and e€Cp(£2), where 
CptfT) denotes the class of piecewise continuously differentiable functions 
in £2". 

We shall denote the boundary of £2 by y. Let 

n ■ l x ' r: nj < °1> 

where |-H. is the outward normal derivative of u along y, y^ is the 
inflow portion of the boundary y. 

Suppose that u satisfies the condition 
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(2.2) inf [max{ |Vu(P | ^ , Au(P)} ] = a > 0, 


p ft 




where j 7u( P ) | 2 = £ 

i=l , n 

(*“ ) 2 . Then 

1 

ft 

can be divided into compact 

subregions ft^ and ft 2 

such that 




| V u | 2 > a 


in , 


Au > a in 


and ft = ft. 

We introduce some notation we shall need later. Let 


u = max u(P), 
max 

P ft 


“rain ’ mIn “ <P> ’ 
P ft 


[u] 


u 

max 


u 


min’ 


s = max 
Pey 



9 


e 


-Au(P) 


PCft , | Vu(P) | 


max { 
“1 


rl- 


We now obtain an a priori bound on the stability of the hyperbolic 
problem (2.1). 
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Theorem 1: Suppose u satisfies condition (2.2). Then for any f for 
which B(u)e = f has a solution e assuming prescribed values along 
the solution is unique, and 

(2.3) Hell < A (u)llfll + D (u)(J |e| p ds) 1/p 

P P P P -y 

T 1 

for all P€ [ 1 ,°° ) • Here 


A 

P 


<■> - ^ 


1/P( P ~ 
P 



crnrr > 


and 


V u> - <srfsy> 1/p ' 

P 


C 0 (u) is defined as 


Cp(u) = min{^- , 1 - i} exp(-pg[u]) if 8 > 0, 


C^(u) = min{i- , 1 - exp(-[u]) i_f g _< 0. 


Proof : Let g(u) be a smooth function of u, which we shall specify 
later. Multiplying equation (2.1) by g(u) | e | p-1 sgn(e) and integrating 
over ft, we obtain 

/ g(u)|e| p 1 sgn(e)f dx = / g(u)|e| P 1 sgn(e) (7»(eVu))dx. 
ft ft 


(2.4) 
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Here sgn(e) is the function defined as 


sgn(e) = 


-1, e < 0, 

0, e = 0, 

1, e > 0. 


Integrating the r.h.s. of (2.4) by parts gives 


(2.5) 


/ g(u) | e | P 1 sgn(e)(V* (eVu))dx = -/ |Vu| 2 |e| P dx 


-/ (p ~ l)(g(u)Vu)*(je| P 1 sgn(e)Ve)dx + / g(u)Je| P ds, 

n y 


(2.6)-/ (p - l)(g(u)Vu)» ( | e j p 1 sgn(e)Ve)dx = -/ (-2 -) (g(u)Vu)» V( | e| p )dx. 


Integrating the r.h.s. of (2.6) by parts, we obtain 


-/ 0^— — )(g(u)7u)» ( | e | p 1 sgn(e)Ve)dx = (-2 -) J il. |Vu| 2 |e| P dx 


p a du 


(2.7) 


+ ( P / (g(u)Au) j e | P dx - ( P -- r- 1 ) / g(u) |ej P ds. 


p ' J 3n 

Y 


Combining (2.5) and (2.7) gives 
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/ g(u)|e| P 1 sgn(e)(V»(eVu))dx 
ft 


/ [“ | |Vu| 2 + (1 - i)g(u)Au] |e| P dx 

ft F F 


( 2 . 8 ) 



(g(u) |^)|e| P ds. 


We wish to choose g(u) so that 


g(u) > 0, and 


[ “ + (1 " ? )g(u)Au] > °» 


for all x€ft. 


We first treat the case 6 > 0, where 8 = max { — Au — . 

P€ft 1 | Vu(P) | 


Let 


g( u ) = - exp(-pg(u - u min )). 


Then 


-1 dg 
p du 


I Vul 2 + 


(1 ~ — ) g(u)Au = [6 


Vu 


+ ( i - I) £i] exp(-p8(u - u m . n )). 


If Au _> a it is easy to see that 

— — — I Vu 1 2 + (1 - — ) g(u)Au > (1 - — )exp(-pg [u] ) . 
p du ' 1 p — p 

O 

Next, suppose Au < a. Then |Vu| > a. From the definition of 8 we 
have 


Hence we obtain 


|Vu| 2 + (1 - ^)g(u)Au > | exp | 


Clearly 


exp<-pB[u)) < g(u) _ 


Next, we treat the case 0 < 0. Then Au > 0 


8(u) -i Mp(-(» - U min )). 


Then 


rSl 7 "l 2 + (1 -F )s(u)Su - ( -ir- +(1 


It is easy to see that in this case 


F 5u l Vu | 2 + o - l)g(u)Au > min (1,1- 


exp(-[ul 2 < g(u) <1 . 


-p0[u]). 


for all x€ 


-) — )exp(-[u] ) 


■i) exp(-[u] ) . 


Also 
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Substituting g(u) in (2.8) we obtain the following inequality 


(2.9) 


where 


Now 


C (u) / |e| P dx + I J (g(u) —■) | e | P ds 

p n p Y dn 

<. 1/ g(u) | e | P 1 sgn(e)f dxdy| , 

£2 

Cp(u) = min((l - ) , —) exp(-p8[u]) if 8 > 0, 

- min(i , 1 - ^-)exp(-[u] ) if 8 < 0. 


|/ g(u) | e | P 1 sgn(e)f dx| < — J |e| P ^fl dx. 
fi a fi 


But by Holder's inequality 


/ | e | P *|f|dx _< (/ | e | ^ p 1 ^ q dx) ^ q ll f II , 
fi fi P 


where — + — = 1 . 

P q 

Hence 


( 2 . 10 ) 


|J g(u)|e| p 1 sgn(e)f dx| < — II ell p ^ q ll f II . 
o ' - a p p 


Applying Young's inequality to the r.h.s. of (2.10), we obtain 
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c (u) 

P ^ „ f C ✓ _P_ ll „ I, P 


|J g(u) I e I P i sgn(e)f dx| _< — ^ — Hell 


( 2 . 11 ) 


c 2 ■'p/q 

c P (u) q ( i )P i 


n f ii 


Thus 


C (u) n c 

■A; llen P < — 

2 p — ap 


.jfcMVj 


(/ |e| p ds) 


Pa^/P 


Ilf II 1 


Since p > 1, this implies 


(j£V /P m p < <^ l/p <J le| p ds ) 1/p + {i (c-f^-) 1 /q (i) 1/p } '««, 


Hence 


(2 - 12) ,e, P < (drW ) 1/P </ M p <“> 1/p ♦ (? <i> 1/p (i> 1/q c-ot ,f V 


And this gives us the required result. 




Remark 1 : It is interesting that the continuous dependence estimate 

(2.3) breaks down for p = °°. Richter [2] proved that the result holds 
for p = 00 . Combining Richter's results with Theorem 1, we conclude that 

e. depends continuously on f in L p (£2) for all pe[l,°°], assuming the 
value of e is prescribed along the inflow boundary. 
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We would now like to investigate the rather general situation where at 
least one of the coefficients of the hyperbolic problem is nonvanishing at 
each point of 12, i.e., we wish to see . whether the continuous dependence 
results remain valid if we replace condition (2.2) by the condition 


(2.13) inf [max {|Vu(P)| 2 , |Au(P)|}] > 0. 

pen 


First suppose that 


(2.14) inf [max {jVu(P)| 2 , -Au(P)}] > 0. 

P€fi 


Replacing u by -u and f by -f in (2.1) so that it takes the form 


B(-u)e = -f j 

it is trivial to see that Theorem 1 remains valid if the following 
modifications are made: 

(i) inf [max |(Vu(P)| 2 , - Au(P)}] = a > 0. 

pen 


(ii) sup [ - j] = 8 

|v u (p)r 


(iii) 


« 2 = {PC12: - Au(P) > a} 
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(Iv) = {P y: > 0} , i.e., y^ now represents the outflow 

portion of the boundary y. 

Richter notes in [2] that if Au is not always positive or negative 
where Vu = 0 it may or may not be the case that B(u)e = f has a unique 
solution for all f€L°°(fi), with appropriate initial data for e. For the 
sake of completeness we cite two examples from his paper illustrating this. 
In Figures 1 and 2, the curves indicate characteristics of u, with arrows 
pointed in the direction of increasing u. Both configurations have one 
maxima and one minima. 



Figure 1 

The first depicts a situation where ft can be separated along a 

characteristic (any one from y. to y„) into two subregions 

A o 

Jlj and such that 

(2.15) | Vu | or Au > 0 in , |Vu| or - Au > 0 in 
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Thus the corresponding hyperbolic problem can be solved uniquely for e with 
initial data specified along y^ or y^. Clearly the continuous dependence 
result (2.3) remains valid for this situation. 



Figure 2 


This is not so for the second configuration because of the presence of 
characteristics going between the maximum and minimum points. This difficulty 
can be circumvented by cutting the domain across characteristics by a line 

into subregions ft^ and • The resulting e would in general be 
discontinuous at the interface y^. Clearly such a solution depends on the 
choice of y A and hence is not unique. Note that Cauchy data would not be 
required at such an interface. 

3. THE NUMERICAL SCHEME 

Our method uses the results of Theorem 2.1 with p = 2. We describe our 
method for the two dimensional parameter estimation problem restricting 
ourselves to the case where the domain ft is the unit square, i.e., 
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n = [o,i] x [o,i]. 

The results described would carry through for the more general situation 
where the boundary of 0 is piecewise smooth, with some technical 
modifications. 

The two dimensional problem can be formulated as follows: 

Let L be a grid of points 


L = {(x ± , y^): (i,j) €A}, 

where A is a finite index set in 7 % . Given a data set of observations 

{uC^, y )}, {f(x, , y )} for u and f at the points (x, , y.)£L 
J J 1 J 

determine e(x) that satisfies the equation 

(3 *° B(u)e = f^ (e + I7 (e |y> = f > (x,y)c?i, 

assuming that the inflow section of the boundary is void. 

We assume that u, e, and f are smooth and that u satisfies condition 

( 2 . 2 ). 

Divide the square [0,1] x [0,1] into a grid of points 

L h = {< x i>y-j> : (i>j> eA h }» 

where 

A h ~ { ( 1 * 3 ) 1 1 < i < N+l, l<j<N+l}. 

We assume that L, C L. 

n — 
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The grid is depicted in Figure 3. Consider now the case where the inflow 
portion of the boundary ^ $. Let 

Y' = {Pey: |^(P) < S 1 

for any 6 > 0. We assume that we are provided Cauchy data for e at a set 
of points 


B = {(x A , y_. ) : (i,j)CA C A} . 


Here (i,j)€A whenever 6L, y^ey', hut A may contain other elements 
besides. We define A precisely below. 



Figure 3 
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Consider the situation in Figure 3 where y' is the line segment AB 
of the boundary x = 1. Here A and B do not belong to the grid of 

points L^. We extend AB to the smallest possible line segment 
^1 h’ P 1 P 6 case > such that 

(1) nsn.h- 

(ii) the end points of y^ ^ belong to the grid L^. 


Then 

B = {(x.., y_.): (i,j)€A c A} 

where 

A =U {(i,j): (x ,y .)eY, 

h 2 

Here the union is taken over all h such that c L. 

Thus we are provided with Cauchy data for e 

eCxj.y^) = wU^y^), (i,j)€A. 

We assume the function w is smooth. 

Let 


B h = { (x i»yj ): C± , j >€A h } , 


A h = { (i » j ) : 


where 
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In the situation depicted in Figure 3, B h consists of the points Pj 
through Pg . 

From the data set of observations {f„} f° r ^ where 

f ij = f ( x i»yj)» 

we construct the piecewise bilinear interpolant f^(x,y) of f(x,y), i.e., 


fh( x i»7j) ^ij 


and f^ is piecewise bilinear in n. 

We approximate the unknown function e(x,y) by a piecewise bilinear 
function e h (x,y). In case Yj h is not void there would be M grid 

points belonging to B^. Then e^(x,y) is the bilinear spline function that 
assumes the values 

e h ( x i»yj) = w(x ± ,yj) for (i,j)CA h 

and whose values at the other mesh points {(x^,y^): (i,j)€A^\A^} have to 

be determined. 

Let 5^ denote the space of piecewise bilinear functions v^(x,y) where 

V h (x i ,y j ) = V ij f ° r ( 1 «j) e V A h’ and 


v h ( x i»yj) = w(x 1 ,yj) for (i,j)€A h . 
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Let and denote the vectors 

^h = ^ (e ij^ (l,j)6A h \A h » 

— h = (i,j)€A h \A h * 

Clearly there Is a one to one correspondence between the linear spline 
functions e^(x,y), v^Cx^), and the vectors JE^ and _Vj 1 . and Vjj 

o 

belong to a vector space of dimension K = (N+l) z - M. 

Finally, we approximate u(x,y) by its cubic B-spline interpolant 
u h(x,y). The cubic B-spline <J> (x) in one dimension is the function 
sketched below [3]: 

4/6 

<p' - V - 0 

x 




Figure 4 
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Let <J>j(x) = “ j)* Then 


u h ( x » y ) = I X q ± j 


where {q^} are obtained by solving a linear system of equations. 

The piecewise bilinear approximation for the unknown function e(x,y) 
that we choose is the function e^(x,y) that minimizes 


1 1 


(3.2) 


I f i B < u h > v h " f h | 2 dxdy 


0 0 


over all v h eSh* 

Substituting the explicit form of the functions u^(x,y), f^(x,y), and 
v h(x,y) in (3.2), we obtain 

(3.3) /‘ |B(u h )v h - f/dxdy - vJc,V h - 2^ + 

0 0 

where 

C h is a symmetric matrix of dimension K by K depending only on 
u h (x,y); 

G h is a vector of dimension K depending on u h (x,y), f h (x,y), and 

w(x,y); 

is a vector of dimension K that corresponds to the linear spline 
function v^(x,y); 

and A, is a scalar depending on f h (x,y) and w(x,y). 
n 11 
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We shall prove In Lemma 2 that is positive definite for all 

h <_ hp, where hQ is a positive constant. 

The piecewise bilinear function ^(x.y) that minimizes (3.2) 

i \ a that minimi zes 

,J ' fc V a h 


corresponds to the vector 


4 


- {< 


e . . )} 

ij J 


(i 


(3.4) 


$h v h - 


over all vectors Thus E^ is obtained by solving the linear system 

of equations 


(3.5) 


c h5h " -Sh. 


Since is positive definite, for h small enough, there exists a unique 

solution of (3.5)> and hence the minimization problem (3.2) has a unique 

solution. 

Let e^(x,y) be the linear spline function that corresponds to the 

vector hi- We claim that e^(x,y) converges to the true solution e(x,y) 
2 

in the L norm at a linear rate of convergence. To prove this result we 
need first to prove a lemma. 

Lenina 1: Let u(x,y) be a smooth function that satisfies condition 

(2.2) in the unit square fi. Let u^(x,y) be the cubic B-spline 

interpolant of u(x,y). Then for h £ h^, where hg is a positive 
constant, and all functions v(x,y) that are piecewise continuously 

differentiable in J2 } the following inequality holds; 
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(3.6) llvll 2 < 2A 2 (u)llB(u h )vll 2 + 2D 2 (u)(/ | v(x,y) | 2 ds) 1/2 , 

Y l,h 

where A 2 (u) and D 2 (u) are defined in (2.3). 

Proof : Since u^(x,y) is a cubic B-spline interpolant of u(x,y), 

u h (x,y) is twice continuously differentiable. 

By a standard result in approximation theory, the following inequalities 

(3.7a) |u - UjJ £ K^h\ 

(3.7b) |Vu - VujJ £K 2 h^, and 

(3.7c) |Au - Au^| <_ K^h^ 

hold, for all (x,y )eft. Here the constants Kj through depend on 

higher derivatives of u(x,y). 

Since u(x,y) satisfies condition (2.2), there exists subregions 
and & 2 of n such that |Vu| _> a in Au _> a in For 

any € > 0 we can choose hg so small that 

(3.8) IVUjJ^ a - € in , Au^ a -€ in & 2 


for all h < h^. 
Recall that 


8 = max {- 


pen. 


Au(P) 

| Vu(P) 


rl > 
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[u] = max {u(P)} - min |u(P)}, and 


s = max | • 


We can now choose hg so small that (3.8) and 


Au h (P) 

max {- j} > 8 - e , 


P€il l | Vu h (P) | 


(3.9) 


t u h ] > [u] ~ , and 


max |i^(P)| > s - € 
P€y 


hold simultaneously for all h < h^. Now 


Y' - (FSy: < 6) 


for some 6 > 0. Let y 2 - Y,\Yj. Then 


1 , - ( p ^ ; f£> 6 ) 


If we choose hg smaller, if necessary, we can obtain 


h ( p ) 2 J » for a11 Pe ^2’ 


whenever h hg. Clearly the inflow section of the boundary for 




If we choose e small enough (and hQ correspondingly small), we can 
always arrange 


(3.10) 


A 2 ( u h) < 2A 2 (u), and 


D 2 (u h ) < 2D 2 (u) 


to hold for all h £ h^. We then obtain a cruder version of estimate (2.3) 


(3.6) 


i vn 2 <_ 


2A 2 (u) 8 B(u^)vll 2 


+ 2D, 


><«)(/ 

Y 


v(x,y) 


2 ds) 1/Z 


l.h 


replacing = {P€y : ^ h ( p ) < 0} by y^ h 

r.h.s. in (2.3). 


in the second term on the 



We shall use Lemma 1 to prove that the matrix defined in (3.3) is 

positive definite for all h £ h^, where hQ is a positive constant. 

Lem 2: Let be the matrix defined in (3.3). Then is positive 

definite for all h £ h^, where hg is a positive constant . 

Proof : Since the matrix depends only on u^(x,y), we may, without 


loss of generality, choose w(x,y) = 0 


Then by (2.3) 
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By Lemma 1 


(3.11) 


11 V h H 2 — ' 4 ^ A 2 < ‘ u ^^ 2 H B ^ u h^ v h" ^ 


since v^Cx.y) = 0, for all (x,y)€y 1 h - 
By elementary means one can show that 


(3.12) 


, „ 2 s h ,.T,. 

V h”2 - 24 W 


Combining (3.11) and (3.12) we obtain 


(3.13) 


V h C h V h - 777TT2 V hV 

96A 2 (u) 


The lemma is proved. 


We can now prove the main result of this section. 


Theorem 2: Consider the parameter estimation problem 


(3.1) 


„/ x 9 / 9u n ,9 , 9u^ c 

B(u)e . _ j_) + ^ (e ^) - f 


in the unit square Q , and assume that we are provided with Cauchy data 


for e _on_ 
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e( ’ t ' y)i (x,y)er ljh * “ (x ' 5 °- 

Here u, e, w, and f are smooth functions, and we assume that u satisfies 
condition (2.2). Let e^(x,y) be the piecewise bilinear function that 
minimizes 


(3.2) 


1 

I 


0 


l l B(u h )v h 



over all piecewise bilinear functions v h (x,y)€S h . Then e h (x,y) converges 
to the true solution e(x,y) in the norm at a linear rate of 

convergence. 


Proof : Let e^(x,y) be the piecewise bilinear function e Sj, that 

interpolates e(x,y) at the mesh points {(x^, y ^ ) : (i , j )€A^\ A^} . Then 
e®(x,y) interpolates e(x,y) at all the mesh points {(x^y^): (i,j)£A h }. 

Since ej 1 (x,y) minimizes 

i 1 , 

(3.2) / / |B(u )v - f | 2 dxdy 

0 0 

over all piecewise bilinear functions v h (x,y)eS h , the inequality 

(3.14) (J J |B(u,)e, - f | 2 dxdy) 1//2 _< (J / |B(u,)e* - f, 1 2 dxdy ) l ^ 2 
00 nnn o 0 


holds . 
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By a standard result in approximation theory we have the inequalities 


(3.15a) 



2 


< v> 2 . 


and 


(3.15b) 


llVe - Ve*ll 0 < K h. 
h 2 — 5 


Here and Kg are positive constants depending on higher derivatives of 

e. 


Further, since fg is the piecewise bilinear interpolant of f we have 


(3.16) 


If - 


f h"2 < V 


where Kg is a positive constant depending on higher derivatives of f. 

Using these and the earlier results (3.7) from approximation theory, it is 
easy to show that 


(3.17) 


11 a o 

/ / |B(u, )e^ - f | Z dxdy _< 0(h). 

0 0 


By the triangle inequality 


(3.18) 


IB( “h ><e h " e h>'2 < ,B(u h )e h - f h'2 + ,B(u h )e h - f h"2- 


From (3.14) and (3.17) we conclude that 


IIB(u h )(e h " e h )l! 2 - 0(h) ’ 


(3.19) 
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Since both e^(x,y) and ejj(x,y) interpolate e(x,y) at the points 
{(x 1 ,y j ): (i,j)€A h }, it is easy to see that 

(3.20) / |e* - e h | 2 ds = 0. 

Y l,h 


From Lemma 1 we have that 

(3.6) llvll 2 < 2A 2 (u)IIB(u h )vH 2 + 2D 2 (u)[J | v| 2 ds] 1/2 

Y l,h 

for all functions v that are piecewise continuously differentiable in "ft. 
Using (3.19) and (3.20) we conclude that 

(3.21) lle h - e^» 2 _< 2A 2 (u)llB(u h )(e® - e ^ ) II 2 0(h). 

By (3.15a) 


lie* - ell „ < 0(h 2 ). 
n L — 

Using the triangle inequality once more we obtain the result 


(3.22) 


lie - e^ll 2 <_ 0(h). 
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Remark: If we were to interpolate e(x,y) by piecewise Hermite cubics 

or by a B-cubic spline, our method would converge at a quadratic rate of 
convergence. It is easy to modify the method to obtain higher orders of 
convergence. 

Finally we indicate how the numerical method we have proposed can be 
adapted in case u does not satisfy condition (2.2) but the more general 
condition 

(2.13) inf [max{ jVu(P) | ^ , |Au(P)|}] = a > 0. 

P€£l 

Consider the situation depicted in Figure 2 where u has both a maximum 
and minimum. Let us suppose that f is such that a unique solution e to 
(2.1) exists. We do not need to specify any Cauchy data for e in this 
situation. As before, we cut the domain £2 across characteristics into two 
subdomains £2^ and such that u satisfies condition (2.2) in £2^ 

and condition (2.14) in £2 2 * 

Let 


a(x,y) = e(x,y) for (x.yje^, and 
b(x,y) = e(x,y) for (x,y)e£2 2 . 

We can find a solution a of 2.1 in the domain £2^ without any Cauchy 
data. The same holds for b in the domain £2 2 . Clearly if a solution e 
to problem (2.1) exists in £1 then e must be continuous along y^. 
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Hence we must have 


a(P) = b(P) for all Pey A . 


Let a^ and denote piecewise bilinear approximations for a and 

b respectively. Then a^ and b^ are the piecewise bilinear functions 
defined in and respectively which minimize 

/ I l B <V V h " f h | 2 dxdy + / / |B(u h )t h - f h | 2 dxdy 

°1 ^1 

+ 0 / |t h - v h | 2ds 

Y A 

over all piecewise bilinear functions v^ and t^ defined in 

and respectively. Here 0 is a weighting parameter. 

In case a smooth solution e to (2.1) exists, it is easy to show that 
the function e^ defined as 


1 e h (x,y) = a h (x,y) for (x.y)^, 
e h (x,y) = b h (x,y) for (x,y)C« 2 

converges to the true solution e(x,y) in the norm at a linear rate of 

convergence. 

In a companion paper [4] , we present a multigrid algorithm for parameter 
estimation problems. In this algorithm we seek a finite dimensional 
representation for e which minimizes the l norm of the residuals of a 
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second order difference approximation of (2.1). Extensive numerical 

experiments indicate that the method converges to the true solution in the 
2 

i norm at a quadratic rate of convergence. 
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